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ABSTRACT 

We use Hubble Space Telescope Wide-Eield Camera 3 (iJ^T/WECS) rest-frame optical imaging to 
select a pilot sample of star-forming galaxies in the redshift range z = 2.00 — 2.65 whose multi- 
component morphologies are consistent with expectations for major mergers. We follow up this 
sample of major merger candidates with Keck/NIRSPEC longslit spectroscopy obtained in excellent 
seeing conditions (EWHM ^0.5 arcsec) to obtain Hct-based redshifts of each of the morphological 
components in order to distinguish spectroscopic pairs from false pairs created by projection along 
the line of sight. Of six pair candidates observed, companions (estimated mass ratios 5:1 and 7:1) 
are detected for two galaxies down to a 3cr limiting emission-line flux of ^ 10“^^ erg s“^ cm“^. 

This detection rate is consistent with a ^ 50% false pair fraction at such angular separations (1 — 2 
arcsec), and with recent claims that the star-formation rate (SER) can differ by an order of magnitude 
between the components in such mergers. The two spectroscopic pairs identified have total SER, SER 
surface densities, and stellar masses consistent on average with the overall z ^ 2 star forming galaxy 
population. 

Subject headings: galaxies: fundamental parameters — galaxies: high-redshift — galaxies: structure 


1. INTRODUCTION 


At redshift z ^ 2 — 3 galaxies are growing rapidly 
and build up a large fraction of their present-day 
stellar mass (e.g. iRnddv et al.l l2QQ8f ). As they grow, 
the increased stellar mass is thought to stabilize 
these systems against gravitational instab ilities result¬ 
ing from their large gas fr actions (e.g., iKassin et al.l 
l2Ql4R^n der Wei et al.ll2Q14D . decreas ing their formerly- 


nign 

tall 

gas-pnase yeiocixy aispersions iiLaw et ai.iizuui d. 
iForster Schreiber et al.ll2009l: iNewman et al.ll2013) 

and causing a morphological transform; 
irregular clumpy starbursts (e.g., 

ation from highly- 
Guo et al.l 120121: 

iLaw et al.1 l2012al: Ivan der Wei et al.l 
ences therein) to the modern-day Hubl 

2014 and refer- 
Die sequence (e.g.. 

iPapoyich et al.l 120051: ILaw et al.l 

I2012bl: IConselicdl2014D. 


One mechanism by which such growth occurs is the 
conversion of massive gas reservoirs into stars. Such 
star formation is observed to occur at a typical rate 
^ SOMc:) yr~^ f o r rest-UV selected galaxy samples (e.g., 
lErb et aI1l2QQ6l: IWuvts et al.l 1201 ill , although this SER 


may represent only a s mall fraction of the gas contin¬ 
ually cycling into (e.g., iGenzel et al.l I2QQ8I: iDekel et al.l 
l2QQ9f ) and out of the galaxies (e.g.. ISteidel et al.ll2QlQf ) 
through large-scale gas flows. Likewise, galaxies also 
grow through both major (mass ratio 3 : 1 or lower) 
and minor (mass ratio 4 : 1 or higher) mergers with 
other galaxies. Such events typically contribute both 
stars and gas, thereby building up the galactic stellar 
spheroid population and providing fuel for future gener¬ 
ations of star formation. The role of mergers and merger- 
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induced star formation compared to in-situ star forma¬ 
tion in building up the present day galaxy population 
has been the subject of considerable deba te, with various 
studies claiming both that mergers ye (Ide Ravel et al.l 
I2QQ9I: [Riiech et al1l2Ql4 iTasca et aI] l2Ql4)~and are not 
(iShapiro et al.l l2QQ8l: IWilliams et all I2Q11I: IWuvts et al.l 
I2Q11I: iKaviraj et al.l l2Q13[ ) maior drivers of star forma¬ 
tion and galactic stellar mass assembly since ^ ^ 4. 

Significant effort has therefore been invested both in 
constraining the evolution of the merger fraction for 
star-forming galaxies (e.g.. iConsehce et al.l I2QQ8L I2Q11I: 
iLotz et aT]l2QQ8L [20111: iRawat et al.ll2QQ8[ ) and in assess¬ 
ing the physical effects of such merge rs on the star for¬ 
mation properties of the galaxies iLaw et ^l2QQ7al 
l2Q12at iLotz et ^ I2QQ8I: iLee et ahl l2Q13f ). One method 
employed by such studies is to use high-resolution 
imaging to quantify disturbances and irregularities in 
the surface brightness profile using a variety of non¬ 
par ametric indices (e.g ., iConselice et al .1120001: ILotz et ^ 
l2QQ4[Law et al.ll2QQ7a[ ). However, it is often challenging 
to intepret such indices unambiguously because 2 : ^ 2 — 3 
galaxies are intrinsically clumpy and irregular and sim¬ 
ilar disturbed morphologies can arise both in merging 
systems and in isolated star form ing galaxies due to inter¬ 
nal dynamical instabiliti es (e.g.. lBournaud fc ElmegreenI 
I2QQ9I: IGenzel et al.ll2Qll[ ). 

An alternative way of identifying major mergers is to 
look for close angular pairs (r < 50 kpc, 6 arcsec at 
2 : ^ 2 —3). When the velocity separation between the two 


components in such a 

pair is < 500 km s ^ (see discus- 

sion bviLin et al.ll2004l 

t iLotz et al.ll2008D numerical simu- 

lations suggest that such systems should predominantly 
trace major galaxy-galaxy mergers during their first peri- 
centric passage and before final coalescence (jLotz et al.l 
I 2 OO 8 L I 2 OIOI). Indeed, when merefer rates deriyed from 

such close pairs (e.g.J 

Bundy et al.l 120091: IWilliams et al.l 

2 OIII: Ide Rayel et al.l 1 

120091: iLoDez-Saniuan et al.l 120131: 

Tasca et al.l 120141) are 

; combined with physically moti- 
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vated timescales for interaction the overall agreement on 
merger rates betw een different studies is relatively good 

One major complication faced by such efforts to con¬ 
strain the merger rate (and the physical characteris¬ 
tics of the merging galaxies) however is the incidence 
of false pairs resulting from chance angular alignments 
of galaxies separated by large cosmological distances. 
In the absence of kinematic information, the false mor¬ 
phological pair fraction can be greater than 50% de¬ 
pending on the adopted impact pyameter (see, e.g., 
discu s sion bvlPatton fc Atfieldll2QQ8l: iQnadri fc Willia^ 
I2Q1QI: iChon et al.ll2Q12l: iLaw et al.ll2Q12al) and whether 
or not photometric redshift selection tec hniques (e.g., 
iKartaltepe et al.l I2QQ7I: iBnndv et al.l l2QQ9f ) can be used 
to help trim the list of potential companions. Although 
corrections for the false-pair contribution can also be es¬ 
timated statistically, such statistical corrections do not 
identify which are the true physical pairs, leading to sub¬ 
stantial uncertainty in the derived properties of major 
mergers as a galaxy class. 

In a recent contribution (jLaw et al.l l2QI2al) we used 
HST imaging data to quantify the major merger fraction 
in rest-UV selected {7Z < 25.5j3 star forming galaxies at 
z 2 — 3 using rest-frame optical morphology. Consis¬ 
tent with studies of similar galaxy samples in the litera¬ 
ture (e.g.. lConselice et ani2QIl[) we found that 23f_l% of 
these galaxies in the redshift range 2.0 < z < 2.5 were 
apparent morphological pairs with physical separations 
of < 16 kpc (^ 2 arcsec), and I6lg% were statistically 
likely to be genuine physical pairs. At such small separa¬ 
tions, the pair morphology would not be apparent in the 
ground-based imaging that forms the backbone of our 
star forming galaxy sample. 

Here we follow up t h e morp hological pair sample pre¬ 
sented by iLaw et M (j2QI2a[ ) with rest-frame optical 
spectroscopy to determine which of the apparent pairs 
have similar spectroscopic redshifts suggesting that are 
likely tc ) coalesce within the next ^ 500 Myr (see Ta¬ 
ble 5 of iLotz et al.ll20I0[) . and to identify whether these 
spectroscopically confirmed pairs have physical charac¬ 
teristics that are any different from those of the rest of 
the star forming galaxy population. 

We give an overview of the parent galaxy sample in 0 
and describe the Keck/NIRSPEC spectroscopic follow¬ 
up observations targeting rest-optical nebular emission 
lines in 0 In 0 we present our results for individ¬ 
ual galaxies and discuss the implications of these re¬ 
sults for the physical characteristics of genuine merg¬ 
ing pairs in ^ Throughout our analysis,we adopt 
a standard ACDM cosmology based on the seven-year 
Wilkinson Microwave Anisotropy Probe (WMAP) re¬ 
sults ([Komatsu et al.l[201 ID in which Hq = 70.4 km s“^ 
Mpc“^, flu = 0.272, and Ha = 0.728. 

2. TARGET GALAXY SAMPLE 

Targets were dra wn from the Keck Bary onic Struc¬ 
ture Survey (KBSS; iTrainor fc Steidel |20I2[) for which 
galaxies are originally identified by UnGlZ color selec¬ 
tion down to IZ = 25.5 and spectroscopically confirmed 
to lie at z ^ 2 — 3 using Keck/LRIS rest-UV s pec- 
troscopy ([Adelberger et al.]l2004lSteidel et al.ll2004l) . As 

^ All magnitudes are given in AB units unless otherwise noted. 


discussed by iConrov et al.l (|2008[) , these galaxies are ex¬ 
pected to typically evolve into ^ L* systems by the 
present day. 

Of the few thousand galaxies in KBSS, 306 lie within 
^ I arcminute of the line of sight to bright background 
QSOs and had rest-frame optical imaging data obtained 
using iJ^r/WFCS as a part of Cycle 17 program GO 
11694 (PI: Law). The details of these observat ions 
have been described at length by ILaw et ahl (|20I2a[) . In 
brief, we used the FI60W (Aeff = 15369 A) filter to 
trace rest-frame optical emission from the target galaxies; 
these data reach a depth of 27.9 AB for a 5cr detection 
within a 0.2 arcsec radius aperture and have a typical 
PSF FWHM of 0.18 arcsec (corresponding to 1.5 kpc at 
z-2). 

We s elected galaxies visually classified by ILaw et ^ 
(|20I2al) as ‘Type IT, for which the rest-optical morphol¬ 
ogy consists of two or more distinct nucleated sources 
of comparable (< 10:1) iLieo magnitude and little to 
no evidence for extended low surface brightness features 
connecting the two components (hence reducing the like¬ 
lihood of multiple clumps within a single low surface - 
brightness disk, e.g. JGenzel et al.kOIlUGno et al.ll20I2f) . 
This selection resulted in 56 galaxies, which we further 
restricted to the subset with secure redshifts in the range 
2 : = 2.00 — 2.65 for which Ha emission lies in the K- 
band between strong atmospheric absorption bands. Ad¬ 
ditionally, we restricted our sample to galaxies with pair 
separations greater than 0.5 arcsec (the minimum sep¬ 
aration at which we can distinguish objects in seeing- 
limited follow-up spectroscopy), and less than 2 arcsec 
(the approximate limit beyond which objects can be un¬ 
ambiguously separated in the ground-based photometry 
that forms the backbone of the KBSS sample, and would 
have been classified as two separate galaxies). We also 
rejected targets whose Ha emission was expected to lie 
extremely close to a bright night sky emission line fea¬ 
ture. In total, seven galaxies met our combined selection 
criteria for the 2012 A observing season, of which we were 
able to observe six targets in the available time. 

The morphologies of these six galaxies are shown in 
Figures [1]-01 we label the primary, secondary, and if ap¬ 
propriate tertiary pieces as ‘I’, ‘2’, and ‘3’ respect ively0 
Assuming that the iLieo magnitude of these pieces is a 
proxy for their total stellar mass, we estimate that the 
mass ratio of these mergers ranges from I:I to ^ 8:1 (see 
Table [T]) . 

As described bv ILaw et al.l (|2QI2aD . iLieo magnitudes 
for these galaxies have been combined with ground-based 
optical imaging (and in some cases ground-based J, K 
infrared imaging and/or Spitzer/IRAC photometry) to 
produce stellar population synthesis models of their spec¬ 
tral energy distribution (SED). The photometry used in 
these models encompasses the light from all components 
of a given galaxyQ and iLieo magnitudes have been cor- 

® Note that primary, secondary, and tertiary features are gener¬ 
ally labelled corresponding to their proximity to the centroid of the 
ground-based object detections and do not necessarily correspond 
with Hiqo magnitude. 

^ Although this may bias the SED fit if the two components are 
unrelated, the ground-based imaging that provides the rest-UV 
photometry is neither deep enough nor sufficiently high resolution 
to reliably determine magnitudes for individual components visible 
in the HST imaging data. This limitation likewise precludes us 
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TABLE 1 
Target Galaxies 


Galaxy 

R.A.^ 

(J2000) 

Decl.^ 

(J2000) 

Secondary Distance'^ 
(kpc) 

Tertiary Distance‘s 
(kpc) 

^160 

(Primary) 

Secondary 
Mass Ratio*^ 

Tertiary 
Mass Ratio*^ 

QI2I7-BXII6 

12:19:31.271 

+49:41:21.90 

14 


24.56 

1:1 


QI2I7-MD16 

12:19:28.407 

+49:40:50.15 

9 


23.44 

4:1 


QI623-BX543 

16:25:57.736 

+26:50:09.44 

9 

7 

22.82 

8:1 

5:1 

QI700-MD103 

17:01:00.321 

+64:11:55.42 

13 


22.51 

2:1 


Q2206-BM64 

22:08:52.360 

-19:43:28.27 

10 

13 

24.03 

1:1 

1:1 

Q2343-BX429 

23:46:22.968 

+12:49:05.55 

10 


24.71 

7:1 



^ Coordinates represent the approximate point midway between the components based on the HST/WFC3 imaging data. 
^ Projected distance of secondary from primary source. 

^ Projected distance of tertiary from primary source. 

^ Estimated from ratio of ifieo magnitudes. 


rected for nebular line emission (see details in iLaw et ahl 
l2Q12a[ ). To fit the galaxy SEDs we use C hariot fc Bruzua l 
(CB13) models in combination with a iChabrieJ (j2QQ3l ) 
initial mass function and a constant star formation his¬ 
tory. 


ondary, and (where applicable) tertiary objects using a 
1 arcsec (7 pixel) wide box, and construct a noise spec¬ 
trum for each target by measuring the rms variations in 
a similar box along blank regions of the slit. 

4. RESULTS 


3. SPECTROSCOPIC OBSERVATIONS 

We obtained spectroscopic observations of these six 
galaxies using the K eck/NIRSPEC longslit spectrograph 
(jMcLean et al.lflo^ with the slit rotated to lie along the 
position angle separating the two components visible in 
the i/Sr imaging. Operating in low resolution mode, the 
NIRSPEC slit measures 42 x 0.76 arcsec, with a spectral 
resolution R ^ 1400 as measured from the widths of sky¬ 
lines, and a detector scale measuring 0.143 arcsec pixel”^ 
along the slit and 4.2 A pixel”^ in the dispersion direc¬ 
tion. Integration times were typically 1 hour per target 
(see Table [2]), composed of four 15-minute exposures be¬ 
tween which the target galaxy was dithered along the 
slit. Each target was acquired by blind offsets from a 
nearby reference star using astrometry derived from the 
HST imaging data. 

Eive of the six targets (Table [2j) were observed on 14 
June 2012 in nearly photometric conditions with near- 
IR iU-band seeing ^ 0.5 arcsec EWHM. In addition, one 
target (Q2343-BX429) was included from previous obser¬ 
vations taken in September 2003 under similar observing 
conditions that serendipitously had the long slit aligned 
with the pair separation axis to within 19° (sufficient to 
incorporate both pieces of the system). 

We reduced the s pectroscopic data u sing a hybrid 
scheme described by iKnlas et al.l (|2012[) that includes 
cosmic ray rejection, image rectification, and a two- 
dimensional sky subtraction algorithm (G. Becker, priv. 
comm) to model the bright OH night-sky emission lines. 
The resulting sky-subtracted spectra are stacked, wave¬ 
length calibrated to the heliocentric vacuum rest-frame 
using the OH emission line features, and flux calibrated 
using observations of the AO infrared standard stars 
HD 1160, HD 203856, and HD 18881 (Vega magnitudes 
K = 7.04, 6.84, and 7.14 respectively). Two-dimensional 
reduced spectra (i.e., slit direction along the abscissa and 
spectral direction along the ordinate) for each of the tar¬ 
get galaxies are shown in Eigures [U - IH Except where 
noted otherwise, we extract spectra of the primary, sec- 

from using photometric redshift estimates to identify p robable false 
pairs fe.g.. IKartaltepe et al.ll2007l : IBnndv et al.|[2009r i. 


We show the HST/F160W morphologies and NIR¬ 
SPEC spectra for our six target galaxies in Eigures [T] 
— m along with estimates of the rms noise spectrum for 
each galaxy. The sensitivity of these spectra is a func¬ 
tion of wavelength; using the rms spectra constructed 
from blank sky regions for each galaxy we estimate that 
the 3(7 limit on spatially and spectroscopically unresolved 
emission ranges from 4 x 10”^^ erg s”^ cm”^ at A = 2.1 
/im to ^ 2 X 10”^^ erg s”^ cm”^ at A = 2.4 jam where 
the thermal background of the Keck/NIRSPEC system 
becomes large. We estimate the uncertainty in central 
wavelength and EWHM of individual lines from boot¬ 
strapped Monte Carlo tests in which realizations of the 
noise defined by the RMS spectra have been added to 
the object spectra. We compute the velocity dispersion 
ay of each target by subtracting off the instrumental res¬ 
olution in quadrature from the measured EWHM of the 
Ha emission line, and estimate the spatial effective ra¬ 
dius of nebular line emission along the slit by fitting a 
Gaussian profile, subtracting the 075 seeing in quadra¬ 
ture from the measured Gaussian EWHM, and dividing 
by 2.36 to obtain the effective la radius. 

4.1. Q1217-BX116 

Q1217-BX116 appears to be a relatively isolated dou¬ 
ble system in which both components have similar H- 
band magnitudes (TJieo = 24.56 for the primary compo¬ 
nent, and 24.74 for the secondary component) and are 
separated by 177 in projection (corresponding to 14 kpc 
at z = 2.19). The secondary component has a circu¬ 
larized effective radius = 1.4 kpc, significantly more 
elongated than the primary {r^ < 0.6 kpc)0 Based 
on SED modeling, Q1217-BX116 is young (202 Myr) 
and low-mass (M^ = 2.4 x IO^Mq) with a SER of 
IIM 0 yr”^. Previous re st-UV spectroscopy obtained 
with Keck/LRIS (see, e.g.. iSteidel et al.l[2010[ and refer¬ 
ences therein) serendipitously included both pieces of the 

^ Circularized effective radii are derived using GALFIT 
(jPeng et al.ll2003j and a model of the HST/WFC3 PSF, and con¬ 
verted to kpc assuming that both piece s lie at the systemic redshift 
of the galaxy. As detailed in lLaw et al.1 (|20I2al L radii are likely sys¬ 
tematically underestimated for sources with iJieo > 24 AB. 
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TABLE 2 

Galaxy Emission Properties 


Galaxy 

ZKa 

Ahck 

(A) 

Primary flux 
(10”^^ erg cm”^) 

Primary 

(km 

Secondary (Tertiary) flux^ 
(10”^^ erg cm”^) 

SLITPA^^ 

N 

''exp 

Q1217-BX116 

2.1940 

20967.5 ±0.2 

3.3 ±0.1 

< 91 

< 0.4 

135° 

4 

Q1217-MD16 

2.6154 

23733.6 ±0.8 

13.2 ±0.8 

103 ± 13 

< 1.8 

184° 

4 

Q1623-BX543 

2.5198 

23105.8 ±0.7 

26.8 ±0.9 

180 ±8 

< 1.8 (4.4 ±0.9) 

205° 

3 

Q1700-MD103 

2.3149 

21760.7 ±0.5 

7.5 ±0.2 

102 ± 12 

< 0.5 

95° 

5 

Q2206-BM64 

2.1948 

20972.4 ±0.1 

4.8 ±0.2 

64 ±4 

< 0.4 (< 0.4) 

64° 

5 

Q2343-BX429 

2.1750 

20842.4 ±0.2 

5.0 ±0.2 

66 ±7 

1.0 ±0.2 

24° 

4 


^ Flux limits are 3cr limits for a spectroscopically unresolved feature at the wavelength of Ho; emission for the primary object. Tertiary 
fluxes or flux limits are only given where applicable. 

^ Redshift is from rest-UV absorption line spectrum. 

Degrees East of North 
^ Exposure time is N x 900 seconds. 


Q1217-BX116 



2.06 2.08 2.10 2.12 2.14 

Wavelength ((am) 


Fig. 1.— Top left panel: HST/F160W imaging data for Q1217-BX116 (z = 2.19); panel measures 7x7 arcsec^. The color map has been 
inverted and uses an arcsinh stretch with the black point set to 28.46 AB pixel”^ (23.0 AB arcsec”^). The orientation of the NIRSPEC slit 
is indicated by parallel red lines. Bottom left panel: cutout of Keck/NIRSPEC 2d spectrum showing the corresponding spatial region along 
the slit and ±2000 km s”^ in the wavelength dimension. The colormap has been inverted and uses a linear stretch. Individual components 
are numbered in both panels, and the red circle in each panel indicates the FWHM of the observational PSF. Right-hand panels: spectra 
of the primary and secondary objects along with the associated RMS noise spectrum; grey shaded bars indicate the wavelengths of strong 
night-sky OH emission features. 


galaxy and showed extremely strong Lya emission with 
no clear absorption features. 

The primary component of the pair is detected in Ha 
to high confidence in our longslit NIRSPEC spectroscopy 
with an observed flux of 3.3 ± 0.1 x 10“^^ erg s“^ cm“^. 
As indicated by Figure [T] (lower left panel), there is no 
evidence of a tilt in the Ha emission profile, which has 
a roughly Gaussian profile along the spatial direction 
with a FWHM 0'.'44 (i.e., consistent with an unresolved 


point source with intrinsic effective radius rua <1-2 kpc 
given the ^ 0'.'5 seeing) and is spectrally unresolved with 
(jy < 91 km s“^. Despite the strong detection of the pri¬ 
mary in the middle of a relatively OH-free spectral region 
^ 700 km s“^ wide, there is no evidence of Ha emission 
from the secondary component. If the secondary compo¬ 
nent is not an unrelated interloper along the line of sight 
to Q1217-BX116, we can use the RMS noise spectrum to 
place a 3cr upper limit of 4 x 10“^^ erg s“^ cm“^ on the 
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Q1217-MD16 
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Ha 
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4 
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0 
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Q1623-BX543 


_1: Primary J 

- fHa = 29.8 ±0.9 J 

X lO '^erg/s/cm^ / 

1 z = 2.5198 " 

2: Secondary 

< 1.8 X 10'^ erg/s/cm^ 


' 1 ' ' ' 1 ' ' ' r ' ' 

_ 3: Tertiary 

- = 4.4 ± 0.9 X 10'^ erg/s/cm^ 

z = 2.5211 [ 

- RMS 

—1—1—1—1—1—1—1— 


HST/F160W 


2 1 


N 




I X I I j 


I I 
I I 
I I 


■ ■ 
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r , lT Iri - I I 
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2.34 


Fig. 2.— As FigurefT] but for Q1217-MD16 and Q1623-BX543 {z = 2.62 and 2.52 respectively). The NIRSPEC spectra for Q1623-BX543 
show the tertiary component, visible in the lower left-hand panel as a small spoke extending to the top right of the primary component. 
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Q1700-MD103 



Q2206-BM64 


HST/F160W 



pin 
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P = 4.8 ±0.2 -12 

X 10 'lJ' erg/s/cm^ - i 
0 
-1 


< 0.4 X 10'^^ erg/s/cm^ -\2 

1 

0 
-1 
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- I 3:Tertiary 


t <0.4x1 O '^erg/s/cm^ - 


rvY'ii''''y^ 
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Wavelength (|jm) 


2.14 


Fig. 3.— As Figure (T] but for Q1700-MD103 and Q2206-BM64 {z = 2.31 and 2.19 respectively). The object on the Eastern end of the 
Q1700-MD103 slit is a foreground star. 
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Q2343-BX429 
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Fig. 4.— As Figure[T] but for Q2343-BX429 (z = 2.17). 

Ra emission flux from a spatially and spectroscopically 
unresolved source (i.e., roughly 10 times fainter than the 
primary). 


4.2. Q1217-MD16 

The two morphological components of Q1217-MD16 
differ by a factor of 3.8 in brightness (i^ieo = 23.44 
and 24.89 for the primary/secondary components respec¬ 
tively) and are separated by I'.'l (9 kpc at 2 : = 2.62). The 
primary and secondary components have similar circu¬ 
larized effective radii of = 0.6 and 0.7 kpc respec¬ 
tively. Similar to Q1217-BX116, SED models indicate 
that Q1217-MD16 is relatively young (80 Myr) and low- 
mass (M^ = 7 X IO^Mq) with a SFR of 88 Mq yr“^. 
Previous Keck/LRIS rest-UV spectroscopy obtained at 
three different position angles suggested two components 
to the interstellar absorption line features at 2 ; = 2.609 
and 2 ; = 2.616, along with Lya emission at 2 : = 2.624. 

As indicated by Figure [21 the primary component is 
well-detected in Ra at 2 ; = 2.6154 with flux 13.2 ± 
0.8 X 10“^^ erg s“^ cm“^, consistent with recent 
Keck/MOSFIRE spectroscopy (Steidel et al. in prep) 
which detected emission at a redshift of 2 ; = 2.6168. 
Although the nebular emission is spatially unresolved 
(fhq; <1.2 kpc after accounting for the observational see¬ 
ing) it has a velocity dispersion of ay = 103 ± 13 km s“^. 
Despite the suggestion of a two-component system from 
the rest-UV spectroscopy, no Ra emission is detected at 
the location of the secondary morphological feature to a 
3a limit of 1.8 ± 10“^^ erg s“^ cm“^. 


4.3. Q1623-BX543 

Q1623-BX543 is a complicated system with three dis¬ 
tinct morphological components within ^ 1 arcsecond 
of each other (Figure |2j) and magnitudes i^ieo = 22.82, 
25.12, and 24.48 for the primary, secondary, and tertiary 
components respectively. The circularized effective ra¬ 
dius of the primary component is 0.9 kpc, comparable 
to the secondary and tertiary features. SED models in¬ 
dicate that Q1623-BX543 is young (9 Myr), low-mass 
{M^ = 5 X IO^Mq), and forming stars extremely rapidly 
at a rate of 515 Mq yr“^ (see discussion in 

Previous i7-band observations of Q1623-BX543 with 
adaptive- optics assisted OSIRIS integral-field spec¬ 
troscopy (|Law et al.ll2QQ9[ ) found that components 1 and 
3 are physically associated, with [O ill] emission features 
offset from each other by 125 km s“^ in velocity and 
6.7 kpc (0'.'8) in projected separation. The present NIR- 
SPEC observations indicate that the primary component 
has an effective Ra radius of ^ 1.2 ± 0.3 kpcfl and a 
broad velocity dispersion ay = 180 ± 8 km s~^ with 
no apparent rotation about a preferred kinematic axis. 
Given the lower spectral and spatial resolution of the 
NIRSPEC data, the measured size and kinematics are 
fairly consistent with the OSIRIS [Olll] A5007 obser¬ 
vations {R ^ 3400, PSF FWHM ^ 0'.'15) which sug¬ 
gested that this component has an effective radius of 
1.1 ±0.1 kpc and a net velocity dispersion of = 153 ±7 


® Uncertainty is dominated by the uncertainty in the observa¬ 
tional PSF. 
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Although the NIRSPEC longslit observations were not 
intended to measure the properties of the tertiary com¬ 
ponent (the position angle of the slit was chosen to cover 
the primary and secondary components) we nonetheless 
detect the tertiary feature at 0'/5 projected distance along 
the slit. Extracting the spectrum from this location gives 
a flux of 4.4 ± 0.9 x 10“^^ erg s“^ cm“^ which is off¬ 
set from the primary by 120 ± 30 km s“^ , consistent 
with the 125 km s“^ derived from OSIRIS [Olll] obser¬ 
vations. The tertiary feature is too faint to obtain a re¬ 
liable estimate of its velocity dispersion, but our derived 
value {(Ty = 1 10^70 km s ^) is consistent with previous 
OSIRIS [Olll] estimates of cTmean = 60 ± 10 km s“^. 
We note that the dynamical mass ratio (8:1) of these 
two components as derived from the OSIRIS data is con¬ 
sistent with the stellar mass ratio estimated using iLieo 
magnitude (5:1; see Table [T]). 

The secondary morphological component located I'.'l 
(9 kpc) in projection to the northeast of the primary was 
not covered by the previous OSIRIS lEU data. While the 
successful re-detection of the tertiary component gives us 
confidence in the ability of the NIRSPEC observations to 
discern faint features down to small angular separations, 
as indicated by Eigure[2]we detect no Ha emission from 
the secondary component to a 3cr limit of 1.8 x 10“^^ erg 

— 1 —9 

s cm . 

4.4. Q1700-MD103 

Q1700-MD103 is significantly larger than most of the 
other galaxies, with primary (secondary) magnitude 
iLieo = 22.51 (23.33) and circularized effective radius 
Te = 3.0 (2.7) kpc. The secondary morphological compo¬ 
nent is highly elongated with a centroid located 1'.'5 East 
of the primary in projection (13 kpc at z = 2.31). Given 
the large effective radius of the primary, it is unsurpris¬ 
ing that the best-fit SED for the gaiaxy corresponds to 
an oid (1 Gyr) and massive (M* = 5 x IO^^Mq) steiiar 
popuiation with an ongoing SER of 49 Mq yr“^. 

Q1700-MD103 has strong Ha emission, with a fiux 
from the primary component of 7.5 ± 0.2 x 10“^^ erg 
s~^ cm“^. While the galaxy is a good candidate for rota¬ 
tion given its high mass fsee. e.g.^lForster Schreiber et al.l 
120091: [Law et al 120091: [Newman et al.1[2MJ) no rotation 
is evident along the axis traced by the NIRSPEC slit. 
Rather, it has an instrumentally-deconvolved velocity 
dispersion = 102 ± 12 km s“^ and a large Ha profile 
of effective radius 2.5 ± 0.2 kpc that roughly matches its 
rest-optical continuum radius. 

In contrast, there is no trace of Ha emission from the 
secondary component to a 3cr limit of 5 x 10“^^ erg s“^ 
cm“^. In this case, we have the benefit of HST/ACS 
E814W imaging (from program GO-10581, PI: A. E. 
Shapiey) tracing rest-frame rsj 2500 A emission at the 
redshift of the primary component to heip us underst and 
this negative resuit. As shown in [Law et al.[ (j2012ai see 
their Eig. 9), the second component is extremeiy red 
compared to the primary (/ 8 i 4 — = 2.7 vs 1.3) and 

These are estimates of the integrated line width from the 
composite spectrum of the galaxy; the OSIRIS data demonstrate 
that the galaxy has a mean internal velocity dispersion of CTmean = 
139 km s“^ (RMS = 32 km s“^ ), with a velocity shear of 39 ± 4 
km s“^ . 


only barely detected in the E814W imaging. Since the 
primary and secondary objects have such different colors 
it is probable that they either lie at different redshifts 
or that the secondary is significantly dustier/older; both 
explanations are consistent with the lack of detection of 
Ha. 

4.5. Q2206-BM64 

Q2206-BM64 is a morphological triple system with 
magnitudes iLieo = 24.03, 23.67, and 23.89 for the pri¬ 
mary, secondary, and tertiary components respectively, 
primary-secondary separation of I'.'l (9.7 kpc) and a 
primary-tertiary separation of 1'.'5 (12.8 kpc). While the 
primary and secondary components are relatively circu¬ 
lar with effective radii = 0.8 and 1.4 kpc respectively, 
the tertiary component is highly elongated with = 1.3 
kpc. Q2206-BM64 has an SED best fit with a relatively 
old (806 Myr) and massive (M* = 3 x IO^^Mq) stellar 
population and 37 Mq yr“^ of ongoing star formation. 

The unusual configuration of all three clumps in a 
nearly straight line meant that it was possible to obtain 
NIRSPEC spectroscopy of all three components simulta¬ 
neously. Although all three have similar colors based 
on our marginally resolved multi-wavelength ground- 
based imaging, as indicated in Eigure [4] only the pri¬ 
mary component is detected in Ha emission with a flux 
of 4.8 ± 0.2 X 10“^^ erg s“^ cm“^. While the emission 
from the primary is spatially unresolved (effective radius 
< 1.2 kpc), it has a velocity dispersion of cTv = 64 ± 4 
km s“^. 

4.6. Q2343-BX429 

Q2343-BX429 is composed of two clumps of magnitude 
Hieo = 24.71 and 26.83 for the primary and secondary 
pieces respectively, with a projected separation of 1'.'2 
(10 kpc at ^ = 2.17). The primary component has an 
effective radius = 1.5 kpc, while the secondary com¬ 
ponent is spatially unresolved with < 0.6 kpc. Sim¬ 
ilarly to Q1700-MD103 and Q2206-BM64, the SED of 
Q2343-BX429 is most consistent with an old (1.1 Gyr) 
and relatively massive (M* = 2 x IO^^Mq) stellar popu¬ 
lation, with a small ongoing SER of 14 Mq yr“^. 

By coincidence, Q2343-BX429 was observed in July 
2003 using Keck/NIRSPEG with the same instrumental 
setup and similar LC-band seeing (0'.'5) as the rest of our 
galaxy targets, and with a position angle sufficiently close 
to the separation vector between the two morphological 
components that both fell within the spectroscopic slit. 
The original redu ction and analysis of t hese data have 
been described bv [Erb et al.[ (|20Ci4 [200^ . As illustrated 
in Eigure m both the primary and secondary components 
are detected in the longslit spectra. The primary com¬ 
ponent has a measured Ha flux of 5.0T0.2 x 10“^^ erg 
s“^ cm“^ with an effective radius < 1.2 kpc and a ve¬ 
locity dispersion cTv = 66 ± 7 km s“^. Although the 
secondary component is faint (Ha flux 1.0 ± 0.2 x 10“^^ 
erg s“^ cm“^) it is centered on exactly the location ex¬ 
pected from the broadband HST/F160W imaging data 
and has a centroid offset from the primary by 140 ± 30 
km s“^. While this secondary component can be de¬ 
tected, we note that its flux calibration has a systematic 
uncertainty (of roughly a factor of 2) compared to the 
primary component since it lies on the edge of the spec- 
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trograph slit. It is also too faint, and too close to an 
OH skyline to obtain a reliable estimate of the intrinsic 
velocity dispersion. 

5. DISCUSSION 

On the whole, the Ha-derived properties of the tar¬ 
get galaxies match well with expectations based on 
HST/WFC3 imaging and broadband photometry. As in¬ 
dicated by Table [31 the Ha and rest-optical continuum 
radii are consistent with each other to within ^ 20%, in¬ 
dicating that the ongoing star formation is spatially co¬ 
incident with the stellar mass accumulated during past 
star formation episodes. This result is limited by our 
inability to measure true two-dimensional profiles from 
the longslit data, but is consistent with previous find¬ 
ings from studies of z ^ 2 star-forming galaxies using 
adaptive-optics assisted integral field spectros copy (e.g., 
[Law et aI]|2QQ9l : iForster Schreiber et al.l[^Qlll ). 

Similarly, the SFR derived from Ha emission line flux 
and global SED modeling are in reasonable agreement. 
We estimate t he SFR from the Ha lu minosity using the 
relation from iKennicntt et al.l (|1994l ) combined with a 
iChabrierl (|2QQ3[ ) initial n iass function a n d an extinction 
correction b ased on the ICalzetti et al.l (1200^ law (see 
discussion bv lErb et al .11201 3 : 1 Rudd v et al.ll2Q12[ ) and the 
SED-derived E{B — V). As shown in Table El extinction- 
corrected Ha and SED-based star formation estimates 
generally agree to within a factor ^ 2 or better. The 
sole exception (Q1623-BX543) has a formal best-fit stel¬ 
lar population age of less than 50 Myr, which is known 
to significantly inflate estimates of the SER above other 
indicators 

Of the six galaxies observed, two were confirmed to 
have nearby companions (Q1623-BX543, Q2343-BX429; 
mass ratios 5:1 and 7:1 respectively) whose projected 
separations and relative velocities indicate that they are 
likely in the process of merging with the central galaxyQ 
As indicated by Eigure O all six primary objects and 
both confirmed companions lie near the star forming 
galaxy main sequ ence for z = 1.5 — 2.5 galaxies (e.g., 
IWnvts et al .1120 111 ). This pair confirmation rate (1/3) is 
consistent with the expectation that ^ 1/2 of apparent 
pairs with angular separations ^ 1-2 arcsec are phys¬ 
ically unrelated superpositions along the line of sight 
based on the statistical distribution of sour ces in the 
HST/W FC3 imaging fields (see discussion by iLaw et all 
l2Q12a[) . This explanation is perhaps particularly likely 
for Q1700-MD103 and Q2206-BM64, for which the pu¬ 
tative companions fall well below the star formation 
main sequence (lower-right open triangles in Eig. [5]). 
However, our results are also compatible with the hy¬ 
pothesis that the SER in merging pairs can differ by 
a factor of 10 between each of the components, even 
when the components have similar rest-optical contin¬ 
uum magnitu des. Such an explan ation was recently 
put forward bv iSchmidt et al.l (|2Q13f ). who used 3F)-HST 
grism spectroscopy to show that apparent morphological 
pairs tended to have nebular emission-line indicators of 
star formation concentrated in just one component of the 

Although as discussed by IPatton fc Atfieldl (|2008lf and 
ILotz et all (| 201 Il f. even close pairs with similar spectroscopic red- 
shifts are not certain to be merging systems since small redshift 
differences can correspond to large separations along the line of 
sight. 




Fig. 5. — Estimated stellar mass and SFR for the primary (filled 
squares) and secondary/tertiary (filled triangles) objects. Stellar 
masses are estimated from the best-fit stellar population model of 
the integrated light of the system (divided up between components 
according to their relative ifieo flux), SFR are derived from the 
nebular emission line fluxes and upper limits. The error bar in the 
lower right corner of the panel indicates the typical uncertainty of 
a given point (which is dominated by systematics). The solid line 
represents the track of the star form i ng ga laxy main sequence at 
2 ; = 1.5 — 2.5 taken from IWuvts et ahl (120111 1. dashed lines indicate 
SFR a factor of two above and below this track. 

pair, ascribing the difference to different gas content in 
the components of the merger. Here we push the SER 
threshold re quired for detection an order of magnitude 
deeper than iSchmidt et Ml (j2Q13l ). but similarly find an 
absence of evidence for pronounced star formation in the 
merging companions. 

In previous analyses of the z ^ 2 star formi ng 
galaxy population we found ([Law et al.l l2QQ7ai l2Q12alla ) 
that rest-frame UV and optical continuum morphol¬ 
ogy was largely decoupled from other physical prop¬ 
erties, with merger-like galaxies having a distribution 
of stellar masses and SER statistically consistent with 
the non-merger p opulation. However, as discussed by 
[Law et Si] (|2Q12aD apparent morphological pairs tended 
to have marginally higher SER surface density than 
the rest of the z 2 star forming galaxy population 
((Esfr) = 10 Mq yr“^ kpc“^ for mergers vs 4 Mq 
yr“^ kpc“^ for non-mergers, with 0.2% confidence in the 
null hypothesis of being drawn from the same distribu¬ 
tion). Additionally, such pairs were much more likely 
to show Lya in emission than the rest of the galaxy 
sample, with many of the strongest Lya sources (e.g., 
Q1217-BX116) ha ving a double-component morphology 
(iLaw et al.ll2Q12c|). This latter trend was also noted by 
I Cooke et al .1 (120101) . whose study of spectroscopically- 
confirmed Lyman Break Galaxy (LBG) pairs at 2 : ^ 3 
found that pairs were much more likely than non-pairs 
to ex hibit Lya in emission (although c.f. iShibnva et al.l 

[ 20 TI . 

While it is impossible to draw statistically robust con¬ 
clusions from our small pilot sample, it is nonetheless 
informative to consider how the properties of our two 
spectroscopic pairs compare to those of the overall star 
forming galaxy population from which they were drawn 
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TABLE 3 

Galaxy Properties 


Galaxy 

™ a 

(kpc) 

(kpc) 

E(B - V) 

M, 

(1O®M0) 

SFRHa" 
(Mq yr-1) 

SFRsed 

(Mq yr-1) 

SFRcomp'^ 
(Mo yr-1) 

Q1217-BX116 

1.4 

< 1.2 

0.17 

2.4 

9±1 

11 

< 1 

Q1217-MD16 

0.6 

< 1.2 

0.21 

7.0 

62 ±4 

88 

< 9 

Q1623-BX543 

0.9 

1.2 ±0.3 

0.30 

5.0 

158 ±5 

515 

23 ±5 

Q1700-MD103 

3.0 

2.5 ±0.2 

0.28 

50 

34 ± 1 

49 

< 2 

Q2206-BM64 

0.8 

< 1.2 

0.21 

30 

15 ± 1 

37 

< 1 

Q2343-BX429 

1.5 

< 1.2 

0.18 

20 

14 ±1 

14 

3±1 


^ PSF-correc ted circula r ized eff ective radius of primary component derived from HST/WFC3 broadband 
imaging. See lLaw et al.l (I2012al) for discussion of typical uncertainties. 

^ PSF-corrected Gaussian half-light radius of primary component derived from Keck/NIRSPEG Ho profile 
along the spectroscopic slit. 

^ Extinction-corrected estimate for the primary component. Uncertainties represent statistical uncertainty in 
the Ho flux measurement, and do not incorporate systematic uncertainty in the SFR prescription adopted. 
^ Extinction-corrected estimate for the morphological secondary companion (tertiary source for Q1623- 
BX543). 


{{ Mg)) ^ IQ^^Mr:^, (SFR) ^ 30 Mq yr~ ^, s ee F ig. 18 
o ^Law et al.l 1201^ . As detailed in 94.31 and 14.61 above. 
Q1623-BX543 has higher than average SFR and Ssfr, 
and lower than average stellar mass and circularized ef¬ 
fective radius. In contrast Q2343-BX429 has lower than 
average SFR, higher than average stellar mass, and 
and XsFR fairly typical of the z ^ 2 star forming galaxy 
population. Neither galaxy shows hya in emission based 
on rest-frame UV spectra. Although our findings are 
therefore marginal ly inc onsistent with the findings of 
iCooke et al.l (|2Q1Q[) and iLaw et aH (|2Q12cD that merg¬ 
ing pairs seem more likely than non-pairs to show Lyg 
emission, they generally match other findings (iLaw et 3] 
l2QQ7ai l2Q12al: iLee et al]l2Q13l: iSchmidt et al.ll2Q13[) that 
spectroscopic pairs do not have physical properties any 
different on average from those of the z ^ 2 star forming 
galaxy population as a whole. 

Looking forward, we note that it will soon be possible 
to resolve such discrepancies. Using catalogs of close-pair 
candidates extracted from large-area morphological sur¬ 
veys such as CANDELS, high-sensitivity follow-up spec¬ 
troscopy using recently-commissioned mnlti-object NIR 
spectrographs su ch as MOSFIRE (e.g., Kiiek et al.ll2Q14D 
and KMOS fe.g.. IWisnioski et~aLli2Q14D is likely to spec¬ 
troscopically confirm samples of a few tens to hundreds of 
genuine mergers for only a modest investment in observ¬ 


ing time. Such statistically large, representative samples 
will allow us to determine the influence of major merging 
events on the evolution of gas and stellar populations in 
galaxies in the young universe. 
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the W. M. Keck Observatory, which is operated as a 
scientific partnership among the California Institute of 
Technology, the University of California, and NASA, 
and was made possible by the generous financial sup¬ 
port of the W. M. Keck Eoundation. DRL and CCS 
have been supported by grant GO-11694 from the Space 
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Inc., for NASA, under contract NAS 5-26555. CCS has 
been supported by the US National Science Eoundation 
through grants AST-0606912 and AST-0908805. AES 
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his IDL-based NIRSPEG data reduction code. Einally, 
we extend thanks to those of Hawaiian ancestry on whose 
sacred mountain we are privileged to be guests. 
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